Potential Coefficients for Ground Return Circuits 

By W. HOWARD WISE 

This paper is concerned with the effect of the finite conductivity and dielectric 
constant of the earth on the potential coefficient for a 1-wire ground return 
circuit. It has been customary to say that the potential coefficient V IQ is 

p\i = c'2 log p" Ip' , elm units per cm. 

It is generally realized of course that this is just a good approximation to the 
true p\i . To see that it is just an approximation one has only to imagine the 
earth turning into air, in which case the distance p" will eventually cease to 
have significance. The object of this paper is to derive the complete expression 
for pn . 
It turns out that 

p,, = c^[2 log p"/p' + 4(M + iN)] (7) 

where 

g-(ft+z)rV"'cos3'fV'ai 



rt = 4irffu, as in Carson's work on Z12 and in mine on Z12 at high frequencies''* 
fe'" = Vl + i{t - \)/2c\a = s 

t = dielectric constant in electrostatic units 
<r = conductivity in electromagnetic units 

= 10"" to I0~'^ in ordinary soil 
X = wavelength in centimeters 
c = velocity of light, in cm per sec. 

M + iN^ vanishes a.sf—>0,f—' co^e^ «= oro-^ =«. 

Ordinarily 4(A/ 4- i^) will not be an important correction to 2 log p''/p'; but if the 
frequency is high and li or :; is small it can be a worthwhile correction. For 
example, if a .02535 inch wire be thrown out on the ground to be a 2 mc antenna 
and we assume thato- = 10'", e = 15 and A = 3 cm. then, with a for wire radius, 

Pn = c'[2 log 2!i/a + 4(M + iN)] 
= 410.455 + .152 + i .319]. 

\/pii is the capacity to ground. If there were two parallel wires the scalar 
potential at the first wire would be I'l = pnQi + Pi/^i . 

Derivation of the Formula 

WE BEGIN with the wave-function for an exponentially propagated 
current in a straight wire parallel to a flat earth. The wave-function 
for a horizontal current-element dipole has been formulated as an infinite 

'John R. Carson: "Wave Propagation in Overhead Wires with Ground Return," 
Bell Sys Tech. Jour. 5, pp. 539-554, 1926. 

= W. Howard Wise: "Propagation of High-Frequency Currents in Ground Return 
Circuits," Proc. I. R. E. 22, pp. 522-527, 1934. 
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integral by H. von Hoerschelmannl The wave-function for the current 
in the wire is obtained by integrating the wave-functions of the current- 
element dipoles along the wire from minus infinity to plus infinity. It is 
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Fig. 1 

The time factor is e'"'. a, b and c are unit vectors pointing in the x, y 
and z directions. 

R, ^ (x^ + y + {h - zyf 

R, = {x^ + / + (/; + 2)^)^'= 

k - 2x/X 

k\ = eiitt)- — i4iraiua in electromagnetic units 

By supposing e to be measured in electrostatic units we can write 

ki = k-(€ — i2cXff)tt. 

^ H. Von Hoerschelmann, Jahrb. der draht. Teleg. 5, pp. 14-188, 1912. 
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It is assumed that fi is everywhere unity in electromagnetic units. 
I = {v'- k'f''\ m - {v-^ - &'"" 

y = a -\- i^is tlie desired propagation constant 
The electric iield ijaraUel to tlie wire is 

= - io e ^12 - -7- 
dx 

It has previously been shown that- 

■ /,, = 70-12 log p"/p' + MQ ~ m] (3) 

where 

Q — iP = — \ {y/i^ + is- — v)e~"'''' cos y'v dv, 

ts^ Jo * 

Tc,' = w sj a. and y' = y -yj a. 

To get the potential cneflicient for a ground return circuit it is necessary 
to compute the scahir potential. 

r = /.r^(^n.+ |-ny + ^n.) 

\bx dy dz J ^ ^^^ 

= Q Pn 

As in previous work (he propagation constant y is assigned the value ik 
as a first approximation. This is an ideal value for 7 but the following 
considerations make it an imperative choice: (l) to assume that the current 
is propagated down the line with a velocity less than that of light makes the 
integrals very hard to evaluate, (2) to assume that the attenuation is not 
zero on an infinite line amounts to assuming an infinite source of energy and 
makes the integrals diverge. 

It should not be inferred tliat the resulting formulas are necessarily poor 
if the physical system does not closely approximate the ideal one in which 
7 is ik. ik is employed as a convenient first approximation in evaluating the 
correction terms in /12 and P12. Eventually, if there were but one wire, 
one would compute 7 = \/(z + 2ii)(G + i<^/pn), wherein Zn and pn 
have been evaluated with ik for 7, and this would be a second a]jproxima- 
tion to 7. Past experience with tlie second approximation so obtained has 
justified the expectation that it would be a satisfactory final result. Since 
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the integrals diverge if the attenuation is not zero the use of an infinite line 
formula presupposes reasonably efficient transmission. 

Since 11 cc e"'''' we have 

arc k tk tw 

Since — — = — or iki = iuQ or I = cQ this is 
dx dt 

ioik-' I- n. = Qc' [2 log p"/p' + 4((2 - iP)]. (5) 

ax 

We have next to consider 

md _ 2a,(l - t') d r^-^^l r JoMe^"'.-rf. 
fe! az iife^ 32 i„ a* Jo (^ + m)(l + T^) 

^ life ^2 Jc (/ -h w)(/ + T^m) J-« dx 

* 
The infinite integral is 

Jo {I -\- tn){l + T^fn)\_ -M J-» J 

Since Jq{v y/^ -\- 'f) and cos k% are even functions of x and sin Ax is an 
odd function of x 

\ Mv Vx^ + y^)e~'^' dx = 2 j MvVx^ + y-) coskx-dx 

= if V < A 

cos 3- V"^ - A^ -r , / 

= 2 . — if k < V 



and so our integral is 

'** A- (r+m)(; + T^m) / 

or, since P = v^ — k^, 



2ik • """'" "^ ^' 



/■" e"""' cos y - rf/ 

or, if we put I = V Va 
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r^ A / ^ , ■ e ~ l\ -2 



and i(kl — k) = 4x0-0) ( 1 + * - ■ - ) = a^ , 



(' 



2C\(7 ) 

2'k r" g-"''''cos y'v-dv 

Va Jo {v + Vv^ + is-) (f + T^ Vv^ + is^) 

where, as in Q — iP, w' = xv -s/a and y' = y 's/a. 

a 

Noting next tliat t- e^"" '' = ~ ■></ave~"' " 
az 

we have 

jQ d ! f" (1 — T^)e~"'' CO?, y' v-v-dv 

? a^ ^' " "^'^ ^ io (i- + Vi^M^^^) (" + t' Vi^M^) 

= —Qc v-, / — „ , e cosyi"(l — T)vdv. 

Since (1 - Or = v + t^ Vc'' + "- - ^^(V'»'^ + iV + v) 

this is ■ _ 

-ecVa Vv^ + iV - " - - , 2 /^-r ^ N~"'' cos y '-(^»' 

i^ Jq L 1- + r Vv'^ + "J 

5 r f" e-^'" cosy'c-(/c "I 

= e4-4(Q - .-P) + 4 / -^^=p^^,J . (6) 

On adding (6) to (5) we have 

QPn = Qr[2 log p"/p' + 4(M + iN)\ (7) 

where 



f" e"*"'" cos y' c 
Jo V*'- + 15^ + "/r 



M -\- iN vanishes as / — > 0, / — » «; , e — * «> or o- — * « . 
When ^2 — ft is minute the leading terms in the approximation (9) 

for M + iiV are ^ - I - -^ log (p" ^/{kl-k^/2) - i ^. 

An Approximation for M + iN 

It is possible to get series expansions for M '\- iN but those which have 
been obtained do not facihtate computation, A fairly good approximation 
to -M + iN is arrived at as follows. 
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Let ie = u , u = e ^ 
e — i2cKa = a — l/r^ 

g' = (// + s - i^f \/a 

g"= {h + z^iy)jVa 

g = (// + 2)fVa 

Since it' + «'')''' = {t' + 2/« + «^ ^ 2/«)^'^ 

= / + « - /«/(/ + «) + ... 
we put 

l/(V/2 + «- + at) = ^/Q ^ u- hi/{l + n) + a/) 
= (/ + ti)/\(a + 1)/^ + (a + l)/« + «^] 
= (/ + ri + r2)/ (a+ !)(/ + rO (/ + r^) 
where fi = «(1 - V l - 4/(a + l)) /2. 
and r2 = m{1 + a/i - 4/(a + l))/2. 
Then 

3/ + .w ^ -'— r(-^ - ^j-'iL!! 



2 

1 



+ '~'Zu 



[~r,e "''■ li (e^"'^') + r^e^^''^ ]i (g-"''^) 



2(a + l)Cr2 - n) 

- r.e""'' li {e-''"''0 + r.e""'' li (e-''"''^]- (9) 
When y = this reduces to 

(a + l)(f2 - ri) 

Jj t .-1 i!t 

where C = .577215665 and, if s = re'", - tt < ^ < tt. 

— s (=0 (—2;' 

The accompanying charts give M and N for 3' = and e = 15. With 
z = li they give M and iV for pu . The computed points are indicated by 
soHd dots on the chart for M; they were obtained iiy numerical integration. 

The approximation (10) was checked against the values obtained by 
numerical integration at a number of points. The discrepancy in each case 
amounted to less than one per cent for both M and A'^. This approximation 
is a much easier way to evaluate the integral than is numerical integration 
but it is a tedious computation with many chances for error. Conse- 




Fig. 2 — M, computed with e ^ 15 and y = 0. The number associated with a curve is the value of (// + s)l;\/a. 
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Fig. 3 — I^, computed with e = 15 and y = 0. The number associated with a curve is the value of (h + s)l;-\/a. 
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quently it is important to observe that a coarser kind of approximation 
may often be good enough. Thus, taking y to be zero, 

Jo u + at a 

If g is very small one might use 



Jo ti + al Ji (1 + aji 



hog(i + ^V-^ii(o. 



Ordinarily precision is not required in M + iN because i(M + iN) is a 
small term in Pi2 . 



